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[57] ABSTRACT 

Apparatus for demodulation and correlation of a code 
modulated 10 MHz signal is comprised of a sample and 


hold analog-to-digital (A/D) converter synchronized 
by a frequency coherent 40 MHz pulse to obtain four 
evenly spaced samples A1,B1,A2 and B2 of each cycle 
of the signal, and means for adding, or subtracting, each 
sample to, or from, one of four accumulators to form 
the sums: 

51, = 2 (.tl - A2) Ml Sla =^(Al- A2) M2 

I “ r 

52, = 2(Bl - B2)Ml S2 q = 1 (Bi - S2) M2 


where 

Ml = 10 MHz Reference • Receiver Code ZO* 

M2 = 10 MHz Reference • Receiver Code Z90* 

The correlation functions that are used for the range 
measurements are then computed from the following 
equations: 

I = 51/COS a — 52/Sin a 
Q ^ 51^cos a — 52^sin a 

where sin a and cos a are determined at the start of 
each range acquisition by forcing the receiver code for 
both Ml and M2 to be equal to 1 and applying these 
equations; 


sin a 


■SI/ 

\si/ + 52 / 


,cos a 


S2, 

Nsi/ + S2/ 


12 Qaims, 8 Drawing Figures 
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FIG. 3 
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DIGITAL DEMODULATOR-CORRELATOR 
ORIGIN OF THE INVENTION 


51, 


Sir 


Nsi/+ 52,^ N 51o^ + S2g^ 


The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 
(72 Stat. 435; USC 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to a method and apparatus for 
demodulation and correlation of a code modulated sig- 
nal with respect to a reference signal. 

The dem^ulation and correlation of a coded signal 
becomes very crucial in many systems, such as in rang- 
ing systems wherein range is essentially determined 
from the time required for electromagnetic waves to 
make a round trip between a transmitter-receiver and a 
tum-aroimd transponder. Ranging systems are primar- 
ily useful in tracking spacecraft, but may also be useful 
in range tracking of aircraft from a ground station or 
another aircraft. In conjunction with simultaneous op- 
eration of three ground stations, enough accuracy may 
be achieved to implement an aircraft collision avoid- 
ance system. Since all of these systems are merely spe- 
cial applications of a ranging system, or are analogous 
to ranging systems in respect to demodulation and cor- 
relation of coded signals, the principles and preferred 
embodiment of the invention will be described with 
reference to a ranging system of the type described by 
Warren L. Martin and Richard M. Goldstein in U.S. 
Pat. No. 3,659,292. 

SUMMARY OF THE INVENTION 

Demodulation and correlation of a code modulated 
signal received with a reference code signal (receiver 
code) is accomplished with automatic operation and 
improved stability by sampling the received signal and 
converting each sample to digital form. Sampling is 
controlled by a frequency coherent signal that is stable 
at a frequency four times the signal to be demodulated 
and correlated such that four equally spaced samples 
Al, Bl, A2 and B2 are obtained at the same phases of 
the signal during successive cycles. Each sample is 
added or subtracted to or from one of four accumula- 
tors to or from the separate sums. 

51, = - A2)Ml Slo=^(Al -A2)M2 

t '‘t 

52, = 2 (Bl - B2) Ml 52„ = 2 (Bl - Bl) Ml 

t “ f 

where Ml is a reference signal at the frequency of the 
received signal multiplied by the receiver code, and M2 
is the reference signal multiplied by the receiver code 
delayed J cycle. The correlation functions are then 
computed according to the following equations: 

/ — S\i cos a — 52/ sin a 

Q = 51q cos a — 52q sin a 

where sin a and cos a are determined in advance by 
forcing the receiver code in the products of Ml and M2 
to be equal to 1, while forming sums in the accumula- 
tors, and computing 


5 


52, 


S2n 


N 51,2 + 52/ \ SlJ+ Slg^ 


As a further feature of the invention, the multipliers Ml 
jQ and M2 are each multiplied by a squarewave chopper 
signal having a period that is long relative to the period 
of the received signal to foreclose contamination of the 
received signal by leakage from either of the other two 
terms of the multipliers. 

‘5 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a functional block diagram of a spacecraft 
sequential ranging system which may employ the pres- 
ent invention to advantage. 

20 fig. 2 illustrates in a graph the inphase, I, and quad- 
rature, Q, correlation functions of two squarewave 
signals which may be the received range code com- 
pared with a receiver code in the system of FIG. 1 for 
ranging. 

25 FIG. 3 illustrates inphase correlation functions for 
three successive codes in the sequential acquisition 
ranging system of FIG. 1. 

FIG. 4 illustrates the functional technique for com- 
paring the received range code with a receiver code in 
30 the system of FIG. 1. 

FIG. 5 illustrates in a functional block diagram the 
demodulation and correlation of the present invention. 

FIG. 6 is a sample timing diagram useful in under- 
standing the present invention. 

35 FIG. 7 is a block diagram of a preferred embodiment. 

FIG. 8 illustrates the organization of an arithmetic 
unit in the preferred embodiment of FIG. 7. 

DESCRIPTION OF PREFERRED 
^ EMBODIMENTS 

Before describing a preferred embodiment, the the- 
ory or ranging will be described in order to better ap- 
preciate the improvement of the present invention over 
the basic binary coded sequential ranging system of the 
45 prior U.S. Pat. No. 3,659,292 (commonly referred to as 
the pt-type system). Ranging is essentially a matter of 
determining the period of time an electromagnetic wave 
takes in traversing the round-trip distance from a base 
station to a transponder. Commonly called the round- 
50 trip-light time (RTLT), this period is measured through 
the use of a carrier modulated by some repetitious code, 
termed a range code. Transponded, the code is returned 
to the base station on a return carrier. The returned 
code appears phase delayed from that being transmitted 
55 due to the round-trip travel time. This phase delay is the 
basis for the range determination. 

Phase delay translates directly into a measure of dis- 
tance if the code’s period is greater than the RTLT. 
This requirement presents no real problem for short 
60 distances as would be the case in such applications as 
surveying on earth or measuring the distance of aircraft, 
but for measuring the range of spacecraft, where the 
code period is some fraction of the RTLT, there is some 
difficulty because the range measurement is ambiguous. 
65 For example, if the transmitted and returned codes have 
periods of 1 second, and differ in phase by 180°, then the 
RTLT to the spacecraft could be 0.5, 1.5, 2.5 .. . 
(«-f0.5) seconds (for any integer «). This ambiguity 
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results due to repetition of the code during the RTLT to 
a spacecraft. RTLT’s to spacecraft are typically on the 
order of tens of minutes for missions to neighboring 
planets, and hours for missions to more distant bodies 
such as Uranus. Unambiguous resolution of such dis- S 
tances would require codes of comparable periods. 
Both the generation and the phase measurement of such 
codes is impractical. Fortunately, it is possible to make 
an a priori estimate of the spacecraft’s range which 
reduces the distance (phase delay) which must be re- 10 
solved. The code period then need be only greater than 
the uncertainty of the a priori estimate since the number 
of code periods within the predicted distance can be 
computed. Adding the predicted number of code peri- 
ods to the fractional period measured by the ft-type 15 
system results in an unambiguous range determination. 

Although measurement precision is conceptually 
quite distinct from ambiguity, it too is related to the 
code period. However, unlike ambiguity, the precision 
to which the returned code phase can be measured 20 
increases as the code period decreases. For example, if 
a precision of 10 nanoseconds is required and the code’s 
period is 1 minute, then phase must be resolved to one 
part in 10*°. This is a very difficult, if not impossible, 
task. However, if the code period is 2 microseconds, the 25 
same precision requires measurement to within only one 
part in 200. 

The /x-type system balances the conflicting require- 
ments of ambiguity resolution and measurement preci- 
sion through the use of several squarewave range codes. 30 
Related by powers of two, these codes are transmitted 
sequentially starting with the highest frequency, which 
establishes the precision of the measurement, and end- 
ing with the lowest which has minimum ambiguity. The 
preferred embodiment of the present invention makes 35 
available codes from about 8 MHz to approximately 1 
Hz. This gives the system a maximum ambiguity resolu- 
tion of 3 X 10* meters and a precision in the neighbor- 
hood of a fraction of a nanosecond. 

The fi-type system of the preferred embodiment is 40 
operated as part of a deep-space-network (DSN) track- 
ing station. Its relation to other ground based subsys- 
tems and to a spacecraft is shown in FIG. 1. The system 
is comprised of a ground station exciter-transmitter 10, 

X- and S-band receivers 11, 12 and dual channel sequen- 45 
tial ranging circuits 13. The X- and S-band receivers 
also employ doppler extractors 14, 15 to provide four 
doppler rate aiding signals X-DOP (0°, 90°) and S-DOP 
(0°, 90°). However, only the I.F. (10 MHz) signals from 
the X- and S-band receivers are utilized to obtain the 50 
improvement in the range data output of the present 
invention. A spacecraft transponder 16 receives range 
codes on an S-band carrier and returns range codes on 
X- and S-band carriers. 

The ranging process starts with the generation of a 55 
range code in the ranging circuits 13. The range code is 
derived from a 66 or 132 MHz frequency reference by 
successive division by powers of two, and the frequency 
reference is produced by a stable generator 17. The 
range code is phase modulated onto the uplink carrier 60 
via modulator 18 and frequency multiplier 19,.and trans- 
mitted to the spacecraft via a power amplifier 20 and an 
appropriate antenna. The transponder aboard a space- 
craft is phase-locked to the uplink carrier. It multiplies 
the carrier frequency by 240/221 and 880/221 to de- 65 
velop, respectively, S- and X-band downlink carriers. 
Concurrently, the received range code is coherently 
detected, filtered in a 1 MHz pass-band channel and 


hard limited, all in a receiver 22, for demodulating the 
two down-link carriers in X- and S-band transmitters 
23,24. 

The down-link X- and S-band signals are received by 
the receivers 11 and 12 which are phase-locked respec- 
tively to the X- and S-band carriers. These receivers 
provide the 10 MHz I.F. signals modulated with the 
range code to the ranging circuits 13 which indepen- 
dently measure the S- and X-band range. The ranging 
system shown thus has two separate and identical rang- 
ing channels, one for X-band, and one for S-band. Only 
the range measurement for the X-band channel will be 
described. It is to be understood that the description 
applies equally to the S-band channel, except for the 
frequencies involved. 

T^e received range code phase differs from that 
transmitted due to RTLT and frequency change due to 
Earth-spacecraft relative doppler, as well as lesser 
phase and frequency variations due to the transmission 
medium. It is impossible to determine the phase differ- 
ence of two squarewaves whose frequencies are not 
identical. Therefore, before range can be determined, 
the doppler effect must be removed. The preferred 
embodiment uses a second code, termed the receiver 
code, to accomplish this. In all respects this code is 
identical to the transmitted range code. The code is 
based upon the same 66 or 132 MHz reference used for 
the range code. Just prior to completion of a RTLT 
after start of the ranging sequence, the receiver code is 
synchronized to the transmitted range code, i.e., they 
are made identical. Following synchronization, how- 
ever, the receiver code is modified by adding properly 
scaled doppler from the doppler extractor to the 66 or 
132 MHz reference. This process, known as doppler 
rate aiding, causes the receiver code to instantaneously 
become a frequency coherent model of the received 
signal. In other words, it becomes identical to what the 
transmitted range code would be if the transmitted 
range code were, modified by doppler. The necessary 
code phase measurement can now be made at leisure 
since the received range code and the receiver code 
remain in a ftxed phase relationship. The resulting range 
determination is equivalent to the “backward looking” 
time of flight at the instant that the receiver code be- 
comes rate-aided. This instant is called time TO. How- 
ever, doppler rate aiding is not a part of the present 
invention; only a part of the preferred embodiment. 

To mechanize the range code phase measurement, 
two correlation voltages are computed for each chan- 
nel. The first voltage, called an inphase voltage, I, is a 
direct comparison of the received range code with the 
receiver code; The second voltage, called quadrature 
voltage, Q, is an analogous comparison of the received 
code with the receiver code delayed by one quarter of 
a code bit period. FIG. 2 shows how the inphase (I) and 
quadrature (Q) correlation voltages vary as a function 
of T, the phase difference. Given any values of I and Q, 
the following equations may be used to determine t in 
microseconds. 


0 S (j),. < -|- given / > 0, Q > 0 

_ T ■ O - 

“ 4 ^ 1 + 


! 4>c <”■ given / < 0, Q > 0 


( 1 ) 


( 2 ) 
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T 1 
’■ = -r ( 


-continued 
— + 1 ) 




4 ' 0 

Jtt 


given I <0, Q < 0 


(3) 


3ir 

2 


- = f +2) 

S < IT given I > 0, Q < 0 


(4) 


r I 
’■ = T< 


/- (? 


+ 3) 


Where r = 


,N 


16 X 2‘ 
3 X F. 


10 


15 


N = code number (e.g., 1 for 4 MHz to 23 for 1 Hz), 

F5 = exciter S3mthesizer frequency (approximately 44 
MHz) from which the 132 MHz reference is derived by 
multiplying by 3. The ranging circuits first determine 
the phase delay for the highest frequency code to estab- 20 
lish the range measurement precision. The same deter- 
mination on a series of lower frequency codes resolves 
the range ambiguity. Because each code is derived from 
the same binary counter, they are phase coherent. 
Therefore it is not necessary to actuily measure the 25 
phase of any component except the first one. FIG. 3 
shows the inphase correlation functions for three suc- 
cessive range components C], C2 *^*id C3. (The quadra- 
ture function has been omitted for clarity.) Naturally, 
their respective correlation functions are coherent. As- 30 
suming tiiat the original (highest frequency) component 
Cl indicates a phase delay of t, the true (or least ambigu- 
ous) range point can be any of the points marked t,(/ = 

1 . . . R)- ^tting the value t into the range tally, and 
retarding the receiver code by the same amount, moves 35 
the observed point on the correlation function to a peak 
for the highest frequency code (one of the points la- 
beled A, B, C, or D). Because the correlation functions 
are coherently related, the range code Q correlation 
function will now be at either a positive or negative 40 
peak. If the peak is positive, the component can be 
ignored since it is already in phase and makes no contri- 
bution to the range tally. If, on the other hand, it is 
negative then one half of the range code bit period is 
added to the range tally and once again, the receiver 45 
code is delayed by the same amount. This is illustrated 
by assuming the code Ci measurement to be at The 
first retardation would leave the receiver code at A. 
The correlation function of code C2 is negative; there- 
fore the receiver code would be moved to point B after 50 
it is delayed by one half of the bit period of code C2. 
Since the correlation function of the ranging code C3 is 
negative at this p>oint, one half of the bit period of code 
C3 would be added to the tally and again the receiver 
c^e is delayed. This final shift leaves the receiver code 55 
at point D. The code phase delay is then determined to 
be equal to T/{. Note that while the first component must 
be explicitly measured, the succeeding component mea- 
surements consist only of determining the sign of the 
inphase correlation voltage. 60 

So far, only ideal conditions have been discussed. In 
practice, charged particle plasmas, both ionospheric 
and interplanetary, corrupt the received signal. The 
range measurement is affected in two ways due to a 
plasma medium. First the range modulation is decreased 65 
in frequency (thereby corrupting code phase) and sec- 
ond the carrier is increased in frequency by a like 
amount (thereby corrupting doppler). This causes a 


slow drift of the measured phase delay. The measured 
drift, known as differential range versus integrated dop- 
pler (DRVID), can be used to determine the time rate 
of change of the total columnar electron content of the 
ray path. To both measure this drift and redetermine the 
initial phase measurement, the ranging program repeats 
the highest frequency component and remeasures its 
phase delay after the range acquisition is complete. 
Although these remeasurements can be performed at 
the correlation function peak, up to a 5 db improvement 
in performance is obtained by retarding the receive 
code by one eighth of the code period. Note from FIG. 
2 that this is the “equal-power” point where both the 
quadrature and in phase voltages have equal magnitude. 
To further increase the amount of DRVID data, the 
highest frequency component is, at the operator’s op- 
tion, transmitted concurrently with all low frequency 
components, except the second, during the resolution of 
range ambiguity. DRVID data is then available 
throughout the ranging acquisition. 

While the ranging algorithm has remained essentially 
unchanged, the present invention represents a marked 
advance and departure from earlier systems in the range 
code demodulation and correlation circuit. Whereas the 
circuit of earlier systems described in the aforesaid pa- 
tent and other sin^ar systems were purely analog, the 
circuit of the present invention is wholly digital. How- 
ever, the departure is more extensive than that, as will 
be evident from the following general discussion of 
demodulation and correlation techniques, and a descrip- 
tion of the demodulation and correlation circuit of the 
present invention in a preferred embodiment. 

The process by which the returned range code is 
compared to the transmitted code is shown functionally 
in FIG. 4. A 10 MHz I.F. signal is obtained from the 
station receiver. To demodulate the range code, a 10 
MHz reference is adjusted through a variable phase 
control unit 30 to quadrature with the I.F. carrier and 
both signals applied to a mixer 31. The mixer product 
after filtering and amplifying in a filter-amplifier 32, is 
the baseband squarewave code signal. 

Two additional mixers 33, 34 are used to multiply the 
baseband squarewave with the receiver code and the 
receiver code delayed by 90°. If two squarewaves are 
allowed to take on the values +1, and input to the 
mixer, the mixer output takes on the product values ± 1 
according to the following table: 


RECEIVED 

BASEBAND 

CODE 

MIXED 

CODE 

MIXER 

OUTPUT 

-1 

-1 

+ 1 

-1 

+ 1 

-1 

-M 

-1 

-1 

+ 1 

-M 

-1-1 


It is easy to visualize that if the phase difference be- 
tween the two signals is 0°, then the output is always 
1, but if it is 90°, the output is alternately -1- 1 and — 1, 
a squarewave with half the baseband code period. 
Given a phase difference of 180°, the output is always 
— 1. For any other phase, the output is a pulse train 
whose ±1 duty cycle is proportional to the phase dif- 
ference. Integration of the mixer outputs results in val- 
ues related to the duty cycle and integration time. Nor- 
malizing the integrated values by their maximum yields 
the I and Q correlation functions in integrator normal- 
izer 35 and integrator normalizer 36, respectively. Inte- 
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gration also averages the received random noise com- 
ponents in the I and Q functions. Demodulation and 
correlation of the range code is now complete. 

Phase adjustment of the 10 MHz reference is crucial 
to best performance of the system. Any misadjustment 5 
degrades the values of the I and Q functions with re- 
spect to noise. Prior to ranging, this phase control is 
“calibrated” or so adjusted as to null the 10 MHz com- 
ponent in the first mixer output. This is the quadrature 
point required for phase demodulation. TThe corre- W 
spending carrier-phase calibration of the present inven- 
tion is fully automatic, as will be described in the fol- 
lowing discussion of FIG. 5. 

Earlier sequential ranging equipment achieved cross 
correlation of received and local codes at 10 MHz using 
balanced mixers 33 and 34 shown in FIG. 4. Since the 
code had not been demodulated from the I.F. carrier, 
and the system bandwidths were relatively wide, small 
phase instabilities in the analog equipment did not affect 
the group delay. However, the system did require daily 
adjustment of the phase shifter 30 used to align the 10 
MHz reference in quadrature to the receiver’s I.F. car- 
rier, and these adjustments were required to be made 
under strong signal conditions to prevent noise from 
obscuring the desired null. Thus, as the received signal 
became weaker, and the receiver AGC adjusted gain to 
compensate, phase shifts occurred rendering the previ- 
ous phase adjustment improper. While the resultant 
losses were small, they nevertheless represented a sys- ^ 
tern deficiency. Furthermore, phase adjustments were 
but once a day, so changes occurring throughout the 
day would not be detected until the following day. The 
present invention overcomes these and other problems 
as will become evident from the following description 
of FIG. 5. 

Referring now to FIG. 5, the demodulation and cor- 
relation of the range code is digital. Two advantages are 
realized over analog techniques, stability and fully auto- 
matic operation. The important distinction between this ^ 
technique and that represented in FIG. 4 is that the 
incoming I.F. signal is immediately coherently sampled 
and digitized by an A/D converter 40. Discrete samples 
thereby created are processed to produce the same I 
and Q functions previously discussed with reference to 
FIG. 4, but in a distinct manner with greater accuracy, 
even with weaker signals. The samples are routed and 
controlled by 40 MHz, 20 MHz, and 10 MHz signals 
developed by a phase-locked loop (PLL) 41 synchro- 
nized to the 10 MHz reference of the base station. As jq 
shown in FIG. 6, the sample and hold jjeriods of the 
A/D converter 40 is controlled by the 40 MHz signal. 
Four samples, Al, Bl, A2, and B2 are taken for each 
cycle of the 10 MHz I.F. input. Samples Al and A2 are 
routed to path Pi by a switch circuit 42 under control of 55 
the 20 MHz signal. Similarly, Bl and B2 are routed to 
path P2by a switch circuit 43 under control of the com- 
plement 20 MHz of the 20 MHz signal. 

FIG. 6 shows a phase offset, <f>, between the 10 MHz 
station reference and the 10 MHz I.F. input. Sample Al go 
is then taken at the point 4> -1- w/4. If -t- ir/4 is desig- 
nated as a, the following equations may be written for 
each of the samples (ignoring noise): 


A\ = Z> sin [a + M{t)] 

(5) 

.^2 = Z> sin [a + IT + Af(0] 

(6) 

Bl = D sin [a + n/2 + M{t)] 

(7) 


8 

B2 = D sin [a +3ir/2+ Af(l)] (8) 

where: 

D is the signal amplitude in digital form, and M(t) is 
the range code modulation k Sin (wj) where k is 
the modulation index. 

Assuming that the multiplication terms Ml and M2 
applied to multipliers 44 to 47 in FIG. 5 are in each case 
the 10 MHz squarewave signal from the phase-locked 
loop, the effect is to multiply samples Al and Bl by 1 
and samples A2 and B2 by —1. Under noise-free condi- 
tions, the outputs SI/, S2/, Sl^and S2gOf the discrete 
integrators 48 to 51, respectively contain the following 
values over a 10 MHz cycle: 

5i; = SliJ = A\-A2 = 2D sin (o + M(t)) (9) 

S2i = S24 = BI-B2 = 2D cos (a + M(t)) (10) 

where the primed terms designate the unprimed terms 
integrated over one 10 MHz cycle. The outputs, I and 
Q, are computed through transformation blocks 52 and 
53 which are preferably implemented with software and 


are defined by: 


/ = 51/ cos a>52/ sin a 

(11) 

Q s 51 q cos a-52<? sin a 

(12) 

Simple trigonometry yields: 


r = Q = 2DM{t)=^2Dk sin [ra^] 

(13) 


where the primed terms again designate the unprimed 
terms integrated over one 10 MHz cycle. From Equa- 
tion (13) it should be apparent that the range code has 
been demodulated. 

The remaining task is to correlate the demodulated 
range code against the output of the receiver coder. 
First allow the I' Value given by Equation (13) to be 
integrated for time t sufficiently large in integrator 52, 
as follows: 

I — ^ r = 2 2 • Z? • A: • sin [wj-/] s 0 (14) 

For t sufficiently large, I = Q s 0. The next step is for 
the right side of Equation (14) to be multiplied in multi- 
plier 44 by the ±1 squarewave k sin [«r(/-f-T)]. If r = 
0, then I is equal to 

Z2-D-k. 

t 

If, however, t is equal to i of the code period, then 

/ = 2 2 • D • fc sin [ea^] sin [aj + 

= 2 2 • • fc • sin [2<aJ] a 0 for large t 

If T is equal to i of the code period, 

/ = 2 2 • D ■ sin [aj] sin [aj + ir] (16) 

= 2 2 • D • sin [aj] (— 1) sin [ra^] 

= —2 2 • D . fc for large L 
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Observe that if the multiplying squarewave is the re- 
ceiver code, then the values of the correlation curve I in 
FIG. 2 are obtained. The values of the correlation curve 
Q in FIG. 2 can be obtained in accordance with the 
present invention in a strictly analogous way in block 5 
53. 

To mechanize Equation (16), Ml is defined to be the 
10 MHz squarewave reference times the receiver code. 
Similarly M2 becomes the reference multiplied by the 
receiver code delayed by i of the code period. These 10 
products Ml and M2 are developed in blocks 54 and 55 
which are mechanized very simply by exclusive-OR 
circuits. Equation (11) applied to Sl/and S2/yields the 
inphase correlation values I while Equation (12) applied 
to Slgand S2g yields the quadrature correlation values 
Q. Both demodulation and correlation are now com- 
plete. However, one thing remains to be discussed: the 
determination of cos a and sin a for use in Equations 
(11) and (12). First, from Equations (9) and (10), if Ml 
and M2 are the 10 MHz reference: 

51, = 51p = £ 2D sin (a+Af(f)) = 2D sin a, for large t. 

52, = S2g = £ 2D cos (o+X/(/)) = 2D cos a, for large t 

From these equations 
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The determination of sin a and cos a in blocks 56 and 57 
is equivalent to “calibration” of the station reference 
previously done by adjustment of the 10 MHz reference 
to a null of the 10 MHz component in the first mixer 
output, but determination of sin a and cos a is done 
automatically under computer control at the start of 
each range acquisition by forcing the receiver code in 
multipliers Ml and M2 to be “1” via switches 58 and 59. 
The resulting SI and S2 values are input to the blocks 56 
and 57 where Equations (19) and (20) are mechanized, 
preferably by software. After “calibration”, multipliers 
Ml and M2 are again selected to be the product of the 
receiver code and the 10 MHz reference. 

It should be noted that noise-free signals were as- 
sumed throughout this discussion. Noise is, of course, 
part of any real signal. Since the sampling of the 10 
MHz I.F. input is sufficiently fast to characterize the 
additive noise in this signal as well as the range-code 50 
signal, the digitally computed I and Q functions become 
a good representation of the inphase and quadrature 
range code correlation values which would be obtained 
from an ideal analog demodulator and correlator. 

An exemplary implementation of the demodulator- 55 
correlator of FIG. 5 for one channel (the S-band chan- 
nel) of the dual channel sequential ranging system 13 of 
FIG. 1 is shown in FIG. 7. The X-band channel is iden- 
tical. 

Since there is only one S-band uplink, only one trans- 60 
mitter coder 60 is required to generate the range code 
directed through a code line driver 61 to the phase 
modulator 18 of the transmitter 10 (FIG. 1). (Tode gen- 
eration is accomplished by dividing the 66 or 132 MHz 
reference from the generator 17 (FIG. 1) by 8 in a di- 65 
vider 62 and applying the resultant 16.5 MHz frequency 
to a 24-stage binary counter in the coder 60. Provision 
is made for a 66 MHz reference to be used to maintain 
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compatability with older DSN equipment. Reference 
will hereafter be made to only the 132 MHz reference to 
illustrate an exemplary embodiment. 

As just noted, the coder 60 consists of 24 flip-flops 
arranged as a binary counter. Repeater flip-flops at the 
output resynchronize the code with 16.5 MHz from the 
divider 62 effectively tying the code to the 132 MHz 
reference, thus promoting stability. 

A receiver coder 63 is preceded by a ^8 pulse adder 
circuit 64. This circuit functions not only to reduce the 
132 MHz to 16.5 MHz, but also to phase shift the re- 
ceiver code. Independent divide by 8 counters provide 
separate control for receiver code relative to the trans- 
mitter code. 

The output of the circuit 64 is used to adjust the phase 
of the receiver code with respect to the transmitter 
code. Adjustment commands may come from either a 
data processor (DP) 65 or from a doppler conditioning 
circuit 66 which receives doppler information from the 
doppler extractor 15. The circuit 64 effectively serves 
to delete a clock pulse in the retard (RTD) mode, or to 
delete a counter state, thus shortening the count se- 
quence, in the advance (ADV) mode. 

Because the 132 MHz reference is derived from the 
44 MHz transmitter frequency synthesizer, the code is a 
subharmonic of the transmitted carrier. This relation- 
ship is a condition precedent to using the doppler ex- 
tractor for rate aiding the receiver coder. However, 
because the link through the spacecraft involves fre- 
quency multiplication (240/221 for S-band; or 880/221 
for X-band), the received doppler must be properly 
scaled before it can be used. Scaling is accomplished by 
applying the inverse multiplier and dividing until the 
doppler frequency has been reduced to one correspond- 
ing to the rate aiding frequency (132 MHz in this exam- 
ple). S-band scaling is accomplished by deleting 19 of 
every 240 doppler cycles providing the ratio 221/240. 
To be useful, the doppler frequency must be further 
reduced to that which would be found at the rate aiding 
frequency. Direction of spacecraft travel is determined 
by the phase relationship between the 0° and 90° dop- 
pler signals. In one case the 90° signal leads the refer- 
ence whereas in the other it lags the reference. 

The codes thus developed by the transmitter coder 60 
and the receiver coder 63 are employed for the trans- 
mission of a range code on the S-band up-link and for 
correlation of the range code received on the S- or 
X-band (the S-band in the example of FIG. 7). At the 
S-band receiver 12 (FIG. 1) there is provided a single 
wideband amplifier with automatic gain control to pro- 
vide a 10 MHz input signal to the A/D converter 40 
(shown in both FIG. 5 and in FIG. 7). 

Comparing FIG. 7 with FIG. 5, it is apparent that, in 
the exemplary implementation, the two channels (0° 
and 90°) are retained, and that code correlation is ac- 
complished using high speed arithmetic units 68 and 70. 
Thereafter, numbers representing the degree of code 
correlation are collected in accumulators 71 to 74, and 
transferred at regular intervals to the data processor 
over a data bus. A coherent multiplier 75 performs all of 
the functions of the phase-locked loop 41 to produce a 
40 MHz sampling command signal and a pair of comple- 
mentary 20 MHz signals for alternating data through 
paired accumulators 71 and 72, and paired accumulators 
73 and 74 corresponding to the respective paired accu- 
mulators 48 and 50, and paired accumulators 49 and 51 
in FIG. 5. The coherent multiplier 75 also performs the 
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analogous functions of the multipliers 54 and 55 in pro- 
ducing signals Ml* and M2* according to the following 
logic functions: 

Ml* = CODE 0‘ ® 10 MHz REF ® CHOPPER, 

M2* = CODE 90‘ ® 10 MHz REF ffi CHOPPER, 

where CHOPPER is a low-frequency squarewave sig- 
nal. Both the coherent multiplier and the receiver coder 
provide inputs to the arithmetic units 68 and 70. Since a 
completely digital demodulation scheme is employed, it 
is necessary to mix (exclusive-OR) the 10 MHz refer- 
ence with the receiver code. Hence the first two terms 
in these logic functions Ml* and M2*. The low-fre- 15 
quency chopper signal is included in these functions to 
foreclose contamination of the received signal by leak- 
age from either the receiver coder or the 10 MHz refer- 
ence. 

Unlike earlier systems which had narrow band filters, 20 
the broadband character of this system might allow 
coupling of both code and/or 10 MHz reference 
through the power system into the AGC amplifier of 
the receiver. The result of such leakage is that the sys- 
tem might measure its own delay rather than the desired 
round trip distance to the spacecraft. The chopper sig- 
nal prevents that. 

TTie chopper signal is generated by the computer in 
accordance with a 10 ms station reference, and has a 
50% duty cycle with a total period of 20 ms. The effect 
is to alternately cause the arithmetic unit to add and 
then subtract. Since the computer is generating the 
chopper signal, it knows the current status (i.e., inverted 
or non-inverted) and can decommutate the samples 35 
which are transferred from the accumulators at 10 ms 
intervals. The result is that any leakage occurring prior 
to the chopper is alternately added then subtracted, and 
hence is cancelled while the desired signal from the 
spacecraft is passed. Measurements indicate that the 40 
leakage is more than 60 db below the received carrier 
with the chopper operating. 

FIG. 8 illustrates an exemplary implementation of the 
arithmetic unit 68. The other arithmetic unit is identical. 

Synchronization and delay flip-flops FF, and FFjare 45 
included to align the coder and chopper sample times 
with those of the A/D converter which is sampling the 
receiver’s I.F. carrier. The 4 bits from the A/D con- 
verter are simultaneously strobed into a hold register 80 
which provides inputs to a 4-bit adder 81. During the 
next sample time, these will either add the new sample 
to or subtract it from previous values accumulated in 
output registers 82 and 83. At the conclusion of the 
add-subtract cycle, the result is transferred to the first of 
the two output registers. 

Two registers are required to implement the alternate 
sampling. Assuming an initial condition wherein the 
registers are cleared, the first sample, Al, is added or 
subtracted into the first output register. Thereafter, the 
second sample, Bl, is likewise transferred to the first 
output register while sample Al is shifted to the second 
register. Note that Al is now available at the second 
input to the adder. Thus, when sample A2 is taken, it 
can be combined with sample Al and stored in the first 65 
register while sample Bl is shifted to the second regis- 
ter. This process continues such that one register con- 
tains 


while the other holds 


n 

2 Ai 


2 Bi 


Since the two output registers contain only 4 bits and 
the samples are collected at a 40 MHz rate, additional 
storage is required. The accumulators 71-74 (FIG. 7 ) 
each consist of 20 stage up-down counters. Wten the 
arithmetic operation results in a carry from the least 
significant 4 bits, an appropriate count command is 
generated and supplied to the proper accumulator. 
Three possible alternatives exist. The accumulator may 
add to, subtract from, or retain its present value. This 
decision depends upon the adders carry output, the sign 
of the present sample, and the status of the adder 81 
(add or subtract) as determined by the flip-flop FF2. 
Two control bits Si and S2are generated by count con- 
trol unit 84 and used to select the count mode of the 
accumulator 71 , for example, according to the follow- 
ing table 


S. 

Sz 

MODE 

0 

0 

Hold (stop count) 

0 

1 

Increnient (count up) 

1 

0 

Decrement (count down) 

1 

1 

Hold (stop count) 


The logic of the count control unit 84 can be expressed 
as 


s, = CARRY 


50 


55 


60 


S 2 = Add . SIGN -H Add • SIGN 

Similar control bits are generated for the accumulator 
73. The control bits for both the accumulators are 
stored for one 40 MHz clock period in the count control 
unit 84. During the subsequent clock cycle they are 
used to retain or change the contents of the accumula- 
tors. 

From the foregoing discussion regarding the arithme- 
tic unit shown in FIG. 8, it is evident that, although 
sampling is at a 40 MHz rate, the collection of 


and 


2 Ai 
i=i 


n 

1 Bi 


proceeds at one half that speed. Reference to the 20 
MHz and 40 MHz waveforms on FIG. 6 will clarify the 
operation. Note that samples are taken during the posi- 
tive portion of the 40 MHz period and held while in the 
negative part. The arithmetic operation takes place 
during this negative period with the result stable and 
ready to be used by the following positive transition. 
Since 



13 


4,112,497 


If n 

2 Ai and 2 Bi 
i=o i=o 

must be separate in order to compute the 10 MHz I.F. 5 
carrier’s phase, it is necessary to provide two accumula- 
tors. The 20 MHz clock can be used to steer the data to 
the proper accumulator as shown in FIG. 7. During the 
negative half cycle data is guided to accumulator A and 
the positive half cycle data is guided to accumulator B. 10 
However, the above explanation is only illustrative of 
the principal of operation. The existence of a register in 
the count control unit 82 (FIG. 8) has the effect of 
delaying accumulation by an additional 40 MHz clock 
period (the clock period of the flip-flops FFi and FF 2 ). 15 
Such a register provides insurance against timing prob- 
lems which can occur when data is transferred from the 
arithmetic unit to the accumulators. 

Returning to the diagram of FIG. 7, it can be seen 
that accumulator B receiv es the 20 MHz clock whereas 20 
the complement, 20 MHz, is directed to accumulator A. 
The data is thus steered to A or B depending upon the 
clock’s phase. Note also that because these clocks are 
derived from a phase-locked loop, wherein a definite 
and certain phase relationship exists between all clocks 25 
and the 10 MHz reference, it is impossible for the sam- 
ples to be directed to the improper accumulator. 

Data is accumulated continuously in the accumula- 
tors. At 10 ms intervals the computer issues a register 
transfer command causing the most significant 16 bits of 30 
the counter to be transferred to a hold register for input 
to the computer. Only 16 rather than all 20 bits are 
transferred since they provide adequate precision and 
facilitate data transfer to a 16-bit computer, such as the 
Digital Equipment Corporation PDP-1 1/20. 35 

Accumulator contents are transferred to the hold 
registers simultaneously for all channels at the conclu- 
sion of a 10 MHz cycle. Referring to FIG. 6, this means 
that contents of Register A is transferred following 
sample A2n, and that of Register B is transferred fol- 40 
lowing sample B2« where the sample number, n, will be 
the same for both transfers. This amounts to simulta- 
neous sampling and forecloses errors due to nonuniform 
sample intervals or samples representing partial cycles. 

Corresponding bits from all hold registers are bussed 45 
in common for transfer to the computer. Individual gate 
controls on the hold registers allow sequential transfer 
of data from a hold register into the computer via the 
data bus. Numerical data accumulated in these hold 
registers represent the degree of correlation between 50 
the received ranging code and the receiver codes. Using 
this information the phase of the received ranging code, 
and hence the range, can be computed in the computer. 

Although particular embodiments of the invention 
have been illustrated and described herein, it is recog- 55 
nized that modifications and variations may readily 
occur to those skilled in the art, and that the basic con- 
cept of the invention may be practiced in still other 
applications. For example, in telemetry of analog data 
from a remote station, such as an oil field, to a central 60 
station, the central station may transmit a reading code 
to a transponder which returns the reading code after a 
delay proportional to the amplitude of an analog signal 
from a transducer. The central station may demodulate 
and correlate the reading code to determine the delay, 65 
and thus determine the value of the analog signal with a 
high degree of accuracy. For more than one transducer, 
the transponder may return a number of reading codes 
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on carriers of different frequencies, one for each trans- 
ducer. Consequently, it is intended that the claims cover 
such modifications and variations. 

What is claimed is: 

1. Apparatus for demodulation of a code modulated 
signal having a predetermined frequency and for corre- 
lation of the demodulated signal with a reference code, 
comprising 

means for sampling said code modulated signal at a 
frequency four times said predetermined frequency 
and for converting each sample to digital form, 
thereby to produce four equally spaced successive 
samples Al, Bl, A2 and B2 in digital form for each 
cycle of said modulated signal, 
four separate means for accumulating multiplication 
products of Ml and M2 with said samples in digital 
form thereby to produce four sums according to 
the following equations: 

51, = 2 (Al - A2) Ml Slg= 2(Al - A2)MZ 

52, = l(Bl - B2)Ml S2 q=^ l(Bl - B2)M2 

where Ml is a reference signal at said predetermined 
frequency multiplied by said reference c^e, and M2 is 
said reference signal multiplied by said reference code 
delayed one quarter of a cycle, 
means for computing inphase and quadrature correla- 
tion functions I and Q, respectively, from said sums 
in accordance with the following equations: 

I = SI/ cos a — S2/ sin a 
Q = Slg cos a — S2q sin a 

where sin a and cos a are determined at the start of 
demodulation operations according to the following 
equations: 


sin a = 


■51/ 

Msi,^ 4- 52/ 


COS a 


52/ 

Nsi/ + S2,^ 


while both multipliers Ml and M2 are made to be equal 
to said reference signal, and 
means for making both multipliers Ml and M2 equal 
to said reference by setting said reference code and 
said reference code delayed one quarter of a cycle 
equal to 1 at the start of demodulation and correla- 
tion operations while sin a and cos a are being 
determined. 

2. Apparatus as defined in claim 1 wherein said means 
for sampling said code modulated signal is comprised of 
means for multiplying said reference signal by four and 
means responsive to the resulting higher frequency 
signal for timing the samples so that said four equally 
spaced samples all occur at the same phases on succes- 
sive cycles of the code modulated signal. 

3. Apparatus as defined in claim 2 wherein said four 
separate means for accumulating samples in digital form 
are synchronized by a timing signal at half the fre- 
quency of said higher frequency signal, and means re- 
sponsive to said timing signal for routing samples Al 
and A2 to both means for accumulating the sums SI/ 
and SIq, and means responsive to the complement of 
said timing signal for routing samples Bl and B2 to both 
means for accumulating the sums S2/and S2g. 
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4. Apparatus as defined in claim 3 wherein said timing 
signal, and its complement, and said higher frequency 
signal, are coherently produced by said frequency mul- 
tiplying means from said reference signal to maintain 
phase coherence of all sampling and accumulating oper- j 
ations with respect to said four samples of each cycle of 
said code modulated signal. 

5. Apparatus for demodulation of a code modulated 
signal having a predetermined frequency, and for corre- 
lation of the demodulated signal with a reference code, 
comprising 

means for sampling said code modulated signal at a 
frequency four times said predetermined frequency 
in response to a sampling control signal, said sam- 
pling means including means for converting each 
sample to digital form, 15 

means responsive to a reference signal at said prede- 
termined frequency for producing said sampling 
control signal so that four equally spaced samples 
Al, Bl, A2 and B2 all occur at the same phases of 
successive cycles of said code modulated signal, 20 
and 

means for accumulating multiplication products of 
Ml and M2 with said samples in digital form in two 
pairs to form the difference of samples Al and A2, 
and the difference of samples Bl and B2 in each of 
two channels by adding or subtracting each succes- ^ 
sive sample in accordance with a first multiplier 
Ml comprised of the product of a reference code 
and a squarewave signal at said predetermined 
frequency, and a second multiplier M2 comprised 
of the product of said reference code delayed one 30 
quarter of a cycle of said squarewave signal to form 
four sums according to the following equations: 


- A2)Ul SIq=^(,AI - A2)m 

S2, = j; (51 - 52) Ml S2n=l (51 - Bl) M2 
t t 


and means for computing inphase and quadrature 
correlation functions I and Q, respectively, from 
said sums in accordance with the following equa- ^ 
tions: 


Q = Sl(? cos a — S2q sin a 

where sin a and cos a are determined at the start of 
demodulation and correlation operations according to 
the following equations: 


while said reference code and said reference code de- 
layed are both set equal to one. 55 

6. Apparatus as defined in claim 5 wherein each of 
said first and second multipliers are further comprised 
of a squarewave chopper signal of low frequency such 
that each chopper cycle spans a number of cycles of 
said code modulated signal. 

7. Apparatus for demodulation and correlation of a 
binary code modulated signal with a reference binary 
code, said signal being at a predetermined frequency 
greater than the pulse rate of said binary code, compris- 

ing 

means for sampling said code modulated signal at 
four equally spaced intervals of each cycle of said 
code modulated signal, and for converting each 
sample to digital form, whereby four samples Al, 


Bl, A2 and B2 in digital form are obtained at the 
same phases of said code modulated signal during 
successive cycles, 

means for accumulating products of Ml and M2 with 
said samples to form four separate sums SI/, S2/, 
Slgand S2g according to the following equations: 


51, = 2(^1 -.42) Ml Slg=X{Al -A2)M2 

52, = 2(51 - 52) Ml 52 q = 2 (51 - 52) M2 


where Ml is the product of a reference signal at the 
same frequency as said code modulated signal and said 
reference code, and M2 is the product of said reference 
signal and of said reference code delayed one quarter of 
a cycle of said reference signal, and 
means for computing inphase and quadrature correla- 
tion functions I and Q in accordance with the fol- 
lowing equations: 

I = SI/ cos a — S2/ sin a 
Q — SI 4 cos a — S 2<3 sin a 

where sin a and cos a are determined in advance by 
forcing said reference code to 1 for both multipliers Ml 
and M2 while producing said sums and computing 


51J + S2] 


51J - 1 - 52j 


8. Apparatus as defined in claim 7 wherein each of 
said multipliers Ml and M2 is a product of a square- 
wave chopper signal in addition to said reference signal 
and said reference code, and said chopper is of a low 
frequency so that each cycle spans a number of cycles 
of said reference code. 

9. Apparatus as defined in claim 7 wherein said means 
for sampling said code modulated signal is comprised of 
means for multiplying said reference signal by four and 
means responsive to the resulting higher frequency 
signal for timing the samples so that said four equally 
spaced samples occur at the same phases on successive 
cycles of the code modulated signal. 

10. Apparatus as defined in claim 9 wherein said 
means for accumulating products is synchronized by a 
timing signal at half the frequency of said higher fre- 
quency signal, and means responsive to said timing 
signal for routing samples Al and A2 to form the sepa- 
rate sums SI/ and Slg, and means responsive to the 
complement of said timing signal for routing samples Bl 
and B2 to form the separate sums S2/and S2^ 

11. Apparatus as defined in claim 10 wherein said 
timing signal, and its complement, and said higher fre- 
quency signal, are coherently produced by said fre- 
quency multiplying means from said reference signal to 
maintain phase coherence of all sampling and accumu- 
lating operations with respect to said four samples of 
each cycle of said code modulated signal. 

12. Apparatus as defined in claim 11 wherein said 
means for accumulating products is comprised of two 
arithmetic units and four accumulators, one arithmetic 
unit for both sums SI/ and S2/ under control of said 
multiplier Ml and one arithmetic unit for both sums 
Slg and S2g under control of said multiplier M2 with 
routing of each sample to be added or subtracted to 
each of two paired sums SI/ and Sljj under control of 
timing signal and to each of two paired sums S2/and 
S2 q under control of the complement of said timing 
signal. 



